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ladnetloe  of  oloeotto  dr^  has  boon  cltod  oa  a  raaoon  for  the 
rolativaly  faatar  aatellng  apaada  than  axpactad  for  aggragataa  oalng 
Stohaa  aattllng  aquation  with  a  particle  packing  aodol  denoity  (Chaae, 
1979;  Hawley,  1982).  However,  the  inability  to  aeaaure  the  actual 


daaalty  of  Individual  partlclea  In  thaae  atudlea  haa  node  It  neceaaary 
to  neka  groaa  aatlnatea  of  one  or  nore  of  the  Stokealan  paranetera, 
uaually  the  aggregate  danalty.  If  aettllng  aggregates  ware  assigned 
too  swell  a  density  hosed  upon  a  given  density  packing  nodel.  It  would 
appear  that  viscous  reduction  had  occurred  In  order  to  provide 
fastar-than-Stokaa  settling  spsads.  Also,  If  aggregates  settled  and 
tiMo  broke  up  at  depth,  the  aaaller  aggregate  coi^onents  observed 
would  appear  to  have  aattled  faster  than  Stokes  equation  allows. 
Again,  viacona  reduction  could  erroneously  be  used  to  explain  this 
aettllng  bdunrlor.  Our  prograw  objectives  were  to  conduct 
aettllng  chart  sr/denalty-gradlwit  ei^rlwents  to  weaanre  aggregate 
particle  denalty  and  vlacoua  drag  Independently. 

A  key  elswant  of  the  project  waa  the  Independent  nsaaurewHit  of 
particle  dsaaity  by  the  uae  of  hi^  denaity  fluida  (iaotonie  widi  die 
aaoaatar)  in  denaity  gradiaat  coluwna.  Aldwui^  the 
k^agr^pkicallyHMMaured  aattling  velocity  and  Stskaa  agnation 
praMctai  wantaorillottito  aggrogate  danaitiea  of  the  <wder  !•<»  g/wl 
Ctiiiiar  to  tko  kigk-otdar  aggragataa  faand  by  Krona  (1978)),  none  of 


in  tiM  pclaani  of 


iiiHtrtta  af  l*A9  Xhia 


/‘Ai 


fluMSd  and  raplaead  ^  tha  aaoalty  of  taa  aarroaadiat  aaala. 


aeiild  r« 


tha  bm^raacy  daa  to  loaar  daaalty  oeeludad  wttsn  of  tin 


erii^aal  a^pratata*  aad  tha  aattliag  atgragata  aevld  haaa  a  aat 
aagatlaa  haoymcy  An  to  tha  clay  (aentaerilloalta)  eoayoBaat. 
■olographic  Isvaatlgatlm  of  tha  bahaalor  of  tha  aggragatas  ao  they 


txaaaraad  a  rtarp  danalty  Intarfaca  has  daaoaatratad 


aspaeta  of 


aggragata  sattllng  baharior  ahldi  ara  of  aora  iaportaaeo  tbw  tha 
objaetivaa  of  tha  initial  pri^na.  Fora  watar  flvdiiag  has  been 
dMonatratad  holographically  to  occur  and  to  have  significant  affects 
an  tha  sattling  bdunrior  of  aggr^ata  particlaa. 


As  a 


iry  of  tha  work  parfomad  under  this  contract  •  this 


report  alll  dascriba  the  baAavior  of  inorganic  aggragatas  eacountaring 
fluids  of  differing  coaposition  than  those  in  uhidi  they  vara  craatad. 
lalagr^liic  tachaiqaas  hare  bean  used  to  show  that  tha  Intarstitial 
eater  of  tinaa  aggragataa  is  not  tightly  bound.  In  fact*  it  la 
aaabangad  rapidly*  and  suggaats  an  initiation  aachaniaa  for  the 
famtioa  of  eartieal  salt  ftogara  Ai  stratiflad  uatar  coluana. 


fha  aattliag  a 
gattaatai  faatwaat 


eonatreetad  for  tl^  paejaet  eith  tha 


eyllagir  to 


tmaatlaa  faiaaa  and 


tiiav  tha  aattliag  fartlelaa  ta  aahiaaa  taaUaal  aattlim  ealaelty* 


r  ahaiwe  afatppal  atlgi  iaal  paMa  far  halagtaphit 

diPPM'’'W;  W  •IM'  VPmHpiBv 
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Thr««  tfpM  of  •sfrogotM  noro  proporod  la  filMT«d»  abiotic 
artificial  aaaaatar.  Thraa  ataadard  clayot  kaolialta* 
aoataorllloalta*  aad  lllita  froa  tha  U.  S,  Gaologlcal  Sarvay  vara  aadi 
allewad  to  aggragata  la  35ppt  artificial  aaavatar  for  about  oao  aoatb. 
Flocculated  aggragataa  oa  tha  order  of  20-300  Mb  vara  gaatly  vlthdrava 
aad  tatzoAiead  both  to  tha  aattliag  ehaabar  aad  the  daaalty  gradleat 
dMobar. 

dlffaraat  daaalty  aadla  vara  uaad  la  the  gradiaata.  la  the 
laltlal  aetup  ruaa  aad  teat*  aerua  daxtraa  vaa  uaed  to  aave  tha 
artraaaly  npeaalva  aatrisaalde  for  latar  ruaa.  KetrlaaBlde-heavy 
vatar  aolutloaa  laotoalc  vlth  35ppt  seavater  vere  created  vith 
daaaltlaa  aa  high  aa  1.40  g/al.  Bvea  at  theae  high  deaaltlea*  the 
aggragotea  coaalataatly  failed  to  becoae  buoyaat*  evea  thoae  vlth 
daaaltlaa  detaralaad  ualag  Stokea  aattliag  equatloo  (vlth  wuiaured 
vatar  vlacoalty)  to  be  laaa  thaa  1.05  g/al. 

III.  iMBAfg 

To  araaiae  tha  beiiavlor  of  tha  aggragataa  aa  they  aattlad  into 
da  tiighor  daaalty  flalda*  tha  hole^caphlc  idccovaloelaatar  vaa 
agplefai.  A  twe-layar*  iaethanwl*  ieotealc  deaeii^  gradleat  vaa 
ftigardd  tmi  artifielal  aaeeafer  (35fft»  l.OMl  g/il»  23*5*€*  0.92 


•neouiit*r«4  th*  int«rfacc»  they  sloMd  donm.  Aftar  •  1«m 

than  a  aacond*  tha  largar  partlelaa  ptmehad  thrM^  dM  latarfa<Ni» 
laavlng  a  atraaaar  of  loifor^aaaity,  ostrudad  pora  watar  in  tbalr 
waka.  Thla  atraaaar  ahom  up  holographically  alnea  tha  dlffaranea  in 
Indax  of  rafxaetlon  hatwaan  tha  pora  watar  atraaaar  and  aadloa  raaulta 
in  Intarfaranca  pattama. 

An  axaapla  of  thla  bahavlor  la  ahoan  In  Flgura  la-a.  Thla  ma 
tha  aacond  axparlaant  parforaad  with  thla  danalty  gradient*  ao  tha 
original  atap-llka  atruetura  had  broken  down  aonawhat  near  tha 
Interface  due  to  the  nixing  Induced  by  tha  aattllng  partlelaa  of  tha 
flrat  axparlnant  (aone  tan  ninutaa  earlier).  Thaaa  unraconatructad 
holographic  Inagaa  ahow  only  the  diffraction  frlngaa.  An  In-focua 
raeonatmetlon  of  flgura  lb  la  ahown  in  tha  Appendix. 

Tha  flrat  vlalbla  particle  to  raauna  aattllng  after  hitting  tha 
Interface  la  laballad  A.  Thla  particle  la  .0256  nlcrona  dlanatar  and 
aatrudaa  a  low  vlacoalty  conduit  in  Ita  wake  that  influancaa  the 
aattllng  of  partlelaa  B  and  C.  In  addition  to  ehamollalng  tha 
aafcaaguant  partlela  aattllng  patha*  thla  conduit  alto  tanda  to  affact 
partiela  norpbology  and  aattllng  orlantatlon.  Tha  reduead  wlaeoalty 
and  danalty  in  tha  conduit  and  tha  aliganant  of  tha  partlela  long  anaa 
with  it»  pamlta  incraaaad  aattllng  apaad  hy  tha  ahhaagnant  partlelaa. 
Thla  Incresscs  tha  llkallhood  that  partlelaa  will  calUia  with  alowar 
■awing  partlelaa  dawnatraan*  and  can  raaalt  ha  alangata  aggragataa 
with  chair  long  ania  arlaatad  wartlcally  (aham  la  Vlinia  »•  A 
dataUad  asplaaatlaa  af  iha  aattllng  hdhawiar  and  Sbathtth  tiilg 
parcielaa  in  fifnea  1  la  pawwidad  la  tha  iiaanaaian. 


i»  m  «  •tHlfit  tetWMHi  that  of  tha  ufipav  and  lowar 

lipMra.  fhp  atvaMMc  l^t  hr  pavtlela  A  attatdhaa  to  idwat  .2  m 
hatffa  tha  «d^|iBal  iMatfaea  (hboat  •  ffaaatara  ot  partlcla  A).  Xf  «a 
Aateiliita  tha  AioMttar  D  of  o  *2  ea  that  eoatataaf  oair 

aaaaitoT,  ita  fllwofit  wouXd  ba  on  the  otdar  of  .0075  a  to  orfar  toat 
tta  aotoaa  aitdi  that  of  tha  ao3»aa  of  rartlcla  A  (.0256  cm  dtoMtar). 

A  Xi^  tar  (X  632.8  an  to  aaqio)  iMiaotog  through  tha  eaatar  of 
a  aoiaaf it  •  eritodar  (a  *  1.3396)  af  .0075  ea  dtoaatax  uould  exit  0.54 


Intsrfercne*  pMk  (bright  fringes)  iriiich  Is  very  close  to  the  ■sxinni 

value  possible.  For  larger  cylinders  with  siallsr  refractive  Indices* 

the  axial  interference  peaks  decrease*  but  never  by  aore  than  50Z. 

For  cylinders  with  d  <  .0025  cn*  the  aaplltude  of  the  Interference 

peak  decreases  rapidly  toward  background  values.  Snail  destructive 

Interference  fringes  (dark  fringes)  will  occur  at  an  angle  of  about 

0.5*  off  axis  and  downstrean  from  the  cylinders  and  will  be  nost 

apparent  for  cylinders  about  twice  (or  Increments  thereof)  the 

diameter  of  our  model  cylinder  (P  *  7*  14*  20*  ...). 

In  auanary*  suxlawmi  fringe  formation  will  occur  on  our  film 

If  p  A  3.9.  For  lower  refractive  index  contrast  (smaller 
c 

I  B^/n^  'll)*  c  larger  cylinder  (presumably  created  by  a  larger 

aggregate)  Is  required  to  produce  similar  constructive  fringes.  For 

1  3330 

fresh  water  over  salt  water  (e.g.*  m  •  *  7359?^*  maximum 

refractive  contrast  approaches  that  of  our  model  cylinder.  Only 
Iceberg  melt  or  river  plumes  over  deep  water  scenarios  (e.g.*  Amason* 
Magdalena)  might  provide  that  much  refractive  contrast.  However*  less 
contrast  from  larger  diameter  cylinders  (larger  particles)  would 
produce  similar  fringes  to  the  ones  we  have  viewed. 

F.  Particle  Dynamics  Mode)  fo**  Aggrenates  with  Occluded  Water 
Flmahlns 

To  understand  the  dynamics  of  porous  aggregates  whldi  exchange 
thalr  occluded  untars  at  some  unknown  rate  with  that  of  the 
aurregndlng  modlum*  un  developed  a  model  to  estimate  the  density  of 
the  occluded  untere  as  a  functlott  of  a  flu^Ung  factor.  To  begin*  uu 
require  aeme  definitiens. 


d 


A  vet  primary  particle  la  defined  as  a  particle  contalnijog  no 
voids  but  has  water  bound  In  the  molecular  matrix.  The  density  of 
such  a  particle,  P^,  Is  defined  as 


M  ♦  M 
s  w 

V  +  V 
s  w 


where  M  and  M  are  the  masses  of  the  solid  material  and  water 
s  w 

respectively,  and  and  are  volumes  of  the  same  respective 
materials.  We  have  determined  In  previous  work  that  the  density  of 
our  samples  of  wet,  primary  montmorillonlte  particles  Is  1.77  gm/ml 
(Gartner  and  Carder,  1977). 

As  the  material  aggregates,  internal  voids  are  created,  partially 
enclosing  occluded  water  (ov).  If  this  ow  Is  tightly  bound,  a  static 
density  results  as  defined  by 


M  +  M  ♦  M 

,  ._5 — * - 2* 

.t.tu 


If  the  ow  Is  exchanged  with  the  surrounding  medium,  then  a  dynamic 
density  results  as 


M,  ♦  ♦  M  • 

0  •  If  Ow 

^dynamic  V,  ♦  ♦  ^ow  * 


M  '  is  an  apparant  aass  of  ow  as  definad  by  the  Stokes  aquation  for  a 

Ow 

sphere  of  diameter  D  aettling  with  velocity  in  a  medium  of  density 
and  viscosity  tu  If  aggregates  or  phytoplankton  contain  the 
occluded  water  (or  cytoplasm)  tightly  enough  that  there  is  no 
significant  flushing  or  exchange  of  pore  waters  with  the  surrounding 
isotonic  medium,  then  it  follows  that 


static  ^dynamic* 


Otherwise 


^static  -  dynamic' 


We  can  describe,  to  the  first  order,  the  effect  of  the  flushing 
of  aggregate  pore  waters  on  the  Stokes  settling  speed  by  developing  a 
relationship  to  describe  p^.  in  terms  of  a  flushing  rate  for  the 
particle.  Suppose  that  some  fraction  f  of  the  ow  is  flushed  each  time 
the  particle  settles  a  distance  of  one  diameter  0.  Then  f/D  F  is 
the  flushing  rate  per  unit  distance  settled.  If  p^(Z) 
and  *  dZ)  arc  the  occluded  water  density  values  at  depths  Z  and 
Z  *  H,  and  the  density  of  medium  p.  increases  linearly  with  depth 


5)  p^(Z)  -  p^(0)  ♦  KZ. 


irfiere  K  ia  a  constant  that  describes  the  change  In  p  with  d^th. 


pwT’' 


k  first  ordsr  spproxlaatioii  of  ths  change  In  In  settling  over 

Ow 

s  short  depth  AZ  •  D  can  be  expressed  as 

^‘*0*^*^  .  **0*^*  ♦  AZ)  -  Pj^(Z) 

AZ  aZ 

fp  (Z  ♦  D)  ♦  (1  -  f)p^(Z)  -  p  (Z) 

St  ^ 

D 

flp.(O)  ♦  K(2  ♦  D))  -  *P«„(2) 

^  _ _ ■ _ _ _ _ ?!? _ 


where  the  fraction  f  of  p^(Z)  that  Is  flushed  Is  replaced  by  the 
denser  water  of  the  asdlun  p^(Z  *  D)  found  at  dapth  Z  D. 
Rearranging  taras*  wu  can  write  the  quatlon  as  a  dlfferantlal  equation 

»««>  *  I  !>«.<*>  ■  I  “  *  *  “> 

or 

Poicz)  ♦  »0«*<Z)  -  m  ♦  r(p^<0)  ♦  KD), 

whara  f  •  f/D. 

TIm  aolatioB  to  this  aqoatlaa  is 


f 


7)  p^(Z)  -  Po^(0)e"”  ♦  ♦  r(Z  +  D  -  l/F)] 


-  P.(0)  -  K(D  -  l/F). 


To  alaulate  our  two-layer,  danaity-gradlant .  particle  dynaalca 
problem,  we  let  P^(0)  *  p^(0)  •  1.024  g/ml  and  the  denalty  eontraet 
gradient.  K  4  .12  g/cn^. 

We  choee  to  exerciae  F  over  a  range  of  valuea  in  order  to 
evaluate  the  effect  that  the  denalty  contraat  of  the  occluded  watera 
with  the  media  haa  on  the  particle  aettllng  apeed  for  varloua 
aggregate  denaltlea.  Aa  mentioned  prevloualy.  the  variation  in 
viacoalty  acroaa  the  interface  had  a  more  elgnlflcant  effect  on 
aettllng  apeed  than  did  the  change  in  denalty.  and  in  order  to  arrive 
at  a  aolutlon.  we  originally  aaaumed  the  denalty  contraat  of  tha 
aggregate  fOj^_  **  p^)  to  be  conatant  with  depth.  Wa  can  mow 
eatimate  how  much  in  error  that  aaaumption  nay  have  been  for  tha 
varloua  particlea  in  thla  atudy. 

How  that  p„  can  be  eatinated  aa  a  function  of  fluahing  rata  F 

’Ow 

and  depth,  one  can  aolve  Stokaa  equation  for  viacoalty  aa  a  faactien 
of  depth,  given  the  aaaaured  aettling  apaeda  of  the  particlea. 


tinea  the  atraamar  from  particle  A  (Figure  1)  affected  the 


viaceaity  r**#  cylinder  in  which  particlea  B  and  C  tinwaled.  it  ia 


important  to  analyse  the  wiaeeeity 


itecad  hy  particle  A  fitet. 


It  ia  claar 


at  .2  cm  helew  tha  erigiaal  iatavfaee  (paattlaa  ef 


particle  A  in  Pigace  la)  thM  a  gradiant  te  rafmctiua 


uatara 


b*lew  th»  iBtcrfaM  m  elM  4«ptli  at  which  the  dcaclty,  vlMMlty,  and 
rafractlwa  ladm  af  tha  Mdiw  with  d^th  (ULaanaaxad 

(■adf—  af  vmraly  Uyar  2  tyya  prapartlaa).  Thua, 

K  •  (1.054  -  i.024)/0.25  cw  -  .12  g/ew^. 

•)  Ml  F)  -  **^^^*^^  “  ^,)Po,(g»y)  -  P,(g)10« 


mi  a^tt) 

•  oM  ff^  la  tha 
if  wUi  aataat.  fa 


5, 


fnai  agaatlaaa  7 
ftaatlaa  af  tha  afpeatati 
ly  dHa  aac  af  agaatiaaa, 
a  walaaity  SMaan 
Is  iOaat  .0745  pafaat 


t  w^CD  la 


alt  » .tl 


t 


of  O.S/cB  and  0.2/ai»  rMpaetlaaly.  TIm  fltuhiag^ilaeter  had  littla  or 
no  affact  on  partielaa  A  and  C  alaca  ■ontaorilloBlta  aada  «p  largar 
fractlona  of  thair  aaaa  valuaa  than  for  partlela  A  (62Z  and  15X  varaoa 
8.6Z»  raapactlvalj. 

Tha  dotted  linaa  rapraaont  ragiona  ohara  two  difforant  aattllag 
tranaltlooa  oeeurrad:  i)  in  Fignra  lb  and  Ic  partlela  A  can  ha  aa«i 
to  ba  aattling  along  or  Juat  adjacent  to  a  atraanar  fr«i  ma  earlier 
particle,  and  11)  In  Figure  Id  and  la,  particle  C  aaana  to  bawa  aovad 
out  of  tha  atraaaar  tuba  of  partielaa  A  and  B.  Particle  A  apparently 
encountered  a  ralatlwaly  eonatant  vlacoalty  whan  In  tha  proxlalty  of 
tha  old  atraanar,  with  tha  vlaeoalty  aharply  Incraaalng  as  It  aovad 
away  froa  tha  atraaaar  raanant.  Particle  C  also  aneountarad  fluid  of 
relatively  constant  viscosity  within  tha  atraaaar  of  particle  A,  with 
aa  abrupt  Increase  occurring  on  paaaaga  out  of  It.  Figure  2  reflects 
thaaa  facts. 

Particle  B  settled  within  tha  atraaaar  of  particle  A  until 
Figara  Id,  where  It  saaas  to  have  edged  out  of  tha  tiAa.  In  Figure  la 
It  appears  to  ba  sll^tly  behind  tha  tuba.  Figure  3  also  corroborates 
this  fact,  slaea  tha  viscosity  aneountarad  by  particle  B  la  clearly 
aach  leas  (within  tha  atraanar  tuba)  than  that  aacountarad  by 
partlela  A  at  aaaa  position  before  tha  atraanar  naa  antmdad. 

Aa  aatirat*  rf  tha  affaetlva  rlacoatty  laecaaaa  aataiia  waraus 
laoida  a  atraanr  is  iadiaacad  bp  iha  anaw  bar  af  ••  caiitii^laa  aa 
Chat  tapraanna  a  rtasaaity  larpnga  'wi 


figara  3. 


■alagiaflqr  of 


oottllog  taroogh  Mooity  iatorfoeoa 


pcooiJoo  •  4ifoct  oboonwtlMi  of  yo»o  wtor  floiiilat  im  to  tlM 
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FlfurM  la  (top)t  lb  («14dla)«  le  (bottoa).  Id  (top:  aaxt  paia),  la 
(■iddia:  aast  paga).  A  aaguaaca  of  holotraaa  takaa  vltb  0.5  aacoed 
apaelag.  Top  arrow  it  idMra  tha  orltlnal  (bafora  axpariwant)  daaalty 
iatarfaca  waa  locatad  aad  rapraaaata  owr  rafaraaca  (Z  •  0)  point. 
Tbraa  aggragataa  labalad  A  (246  m  diaaatar)*  •  (76.5  im  diaaatar). 
and  C  (113.5  ^  dianatar)  ara  tattling  along  tha  aaaa  pathwap*  which 
la  narhad  by  a  bri^t  pora-watar  atraanar  ananating  freai  partlela  A. 
In  rignraa  Id  and  la  partielaa  B  and  C  appear  to  ba  dirartad  fron  thia 
atraanar*  pathapa  balng  inflwancad  by  tha  proxiaity  of  tha  rannant 
atraanar  on  tha  innaillata  laft.  Rota  tha  apparant  fomatlon  of  an 
alangatad  aggragata  in  tha  atraanar  Jnat  above  particle  C. 


Figure  2  (bottoa).  Hologran  of  eevcral  olongatod  aggrogatoa  found 
about  1  ca  below  the  Interface  (note  especially  the  lower  left 
comer).  These  streanllned  aorphologies  suggest  their  fomatlon  to 
bare  occurred  In  a  streaaer.  Note  the  continued  presence  of  low  Index 
atreaaers  even  this  far  below  the  Interface,  suggesting  relatively 
slow  flushing  rates. 


Figure  3.  Viscosity  profile  in  the  density  gradient  traversed  by  and 
Inferred  froa  particles  A,  and  C.  The  sudden  Increase  in  viscosity 
observed  for  particle  C  some  .11  ca  below  the  Interface  occurred  when 
the  particle  left  the  low  viscosity  streamer  and  suddenly  encountered 
aore  resistance  to  settling. 
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In-£ocus  images  of  particles  A,  B,  and  C  after  holographic 
recons  true  t Ion . 
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Image  analysis  techniques  for  holograms  of  dynamic 
oceanic  particles 

Paul  R.  Payne,  Kendall  L.  Carder,  and  Robert  G.  Steward 

A  holographic  image  analysis  system  has  been  developed  to  measure  position,  velocity,  size,  and  shape  of  mi¬ 
croscopic  particles  settling  in  3-D  space.  Images  of  particles  recorded  sequentially  on  individual  holograph¬ 
ic  frames  are  reconstructed  using  an  in-line  far-field  configuration,  registered  in  3-D  space,  and  particle  dis¬ 
placements  (velocities)  between  sequential  frames  are  determined.  Particle  settling  velocity  is  calculated 
using  elapsed  time  between  frames.  Digitized  video  signds  of  the  reconstructed  holographic  images  are  pro¬ 
cessed  to  determine  particle  size,  shape,  and  area  to  facilitate  identification  by  shape  from  frame  to  frame 
and  to  calculate  particle  specific  gravities.  A  cataloging  system  was  developed  to  provide  efficient  data 


management. 


I.  Mrodudfon 

Measurements  of  size,  shape,  and  density  are  critical 
factors  in  the  study  of  particle  dynamics  and  settling 
characteristics  of  marine  hydrosols.  Micrographic 
holography  has  been  widely  in  aerosol  studies*  and 
holographic  movie  cameras  have  been  developed  for 
studying  zooplankton  feeding  behavior  in  ocean  waters.^ 
Using  a  similar  holographic  technique,  the  size  and 
settling  behavior  of  aggregates  resulting  from  mixing 
riverine  water  with  seawater  has  been  determined  in  the 
laboratory^  and  the  effect  of  particle  shape  on  the 
settling  velocities  of  primary  hydrosols  (unaggr^ated) 
has  aim  been  investigated.^-^ 

Carder  et  at.*  have  developed  an  in  situ  holc^aphic 
microvelocimeter  for  the  study  of  microscopic  particles 
suspended  in  seawater.  This  device  records  the  image 
of  each  particle  in  3-D  space  on  a  series  of  successive 
hdograms  with  respect  to  time.  Horizontal  and  vertical 
dimensions  and  cross-sectional  area  have  been  used  to 
determine  particle  density  using  the  Stokes  theorem.^ 
Edge  coordinates  are  useful  for  particle  identification 
(e.g.,  Zemike  moments)  and  analysis  of  shape  on  par¬ 
ticle  rotation  and  settling  speed.^  During  tliM  studies. 
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a  great  number  of  samples  must  be  analyzed  in  order  to 
provide  sufficient  statistics  for  determining  r'-  Jcle 
population  characteristics  of  shape  and  size  as  veil  as 
settling  velocity,  trAjectory,  orientatimi,  and  oscillations. 
Such  measurements  can  be  facilitated  by  use  of  a  com¬ 
puter-controlled  image  analysis  system  to  improve  the 
accuracy  and  reduce  the  number  of  man-hours  in  ana¬ 
lyzing  holographic  particle  data. 

inrowivffi 

A  measurement  technique  was  needed  to  determine 
the  position,  size,  shape,  orientation,  and  velocity  of 
microscopic  partides  moving  in  3-D  space.  Hie  marine 
particles  of  interest  in  the  study^  dtMl  above  typically 
ranged  from  5  to  250  pm  in  dianieter,  whidi  necessitated 
the  use  of  microscopic  techniques  for  any  detailed  size 
or  shape  analysis.  Since  oceanic  particle  densities  may 
vary  from  ~1.03  to  5.2  g/mliter  and  typical  settling  ve- 
lodties  may  range  from  <0.0001  to  1.0  cm/sec,  a  varmty 
of  frame  or  sample  periods  from  <1  sec  to  as  much  as  1 
min  may  be  necessary  in  order  to  provide  short-term 
velodty  measurement  accuracy,  additkm,  the  du¬ 
ration  or  exposure  time  for  the  positicm  measurement 
should  be  <1/500  sec  in  order  to  maintain  suffident 
positional  accuracy  and  prevent  hologram  smearing. 
Therefore,  a  relatively  hi^-speed  mkroacopic  tech¬ 
nique  was  developed*  to  record  the  images  of  multiple 
parties  as  they  settled  in  3-D  HMce.  Todealwiththe 
enormous  data  volvune  genmrated  in  applsring  these 
techniques  (mwe  than  SO  partides/hokigram),  an  au¬ 
tomated  method  for  reconstnicting  and  analyting  the 
partide  data  has  hem  devdoped.  Edge  ooonuiates  are 
recorded  adien  needed  for  partide  fakotiffcatton  and  for 
analgia  of  the  particle  shape  effects  on  partide  rota- 
tim. 
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Fig.  1.  Holographic  image  analysis  system  block  diagram. 


■Ik  vyROTi  uMCfipoon 

A  functional  block  diagram  of  the  holographic  image 
analysis  system  is  shown  in  Fig.  1.  The  major  compo¬ 
nents  and  functions  are:  (1)  a  laser  to  reconstruct  the 
holographic  image,  (2)  a  mechanical  stage  to  position 
the  hologram  in  S-D  space,  (3)  a  magniflcation  lens  and 
rear  projection  screen  on  which  to  reconstruct  the  par¬ 
ticle  image,  (4)  a  TV  camera  to  produce  a  video  signal 

the  image,  (5)  a  video  digitizer  to  convert  the  video 
signal  into  dijptd  data,  (6)  a  TV  monitm  to  display  the 
particle  image,  (7)  a  microcomputer  with  disk  storage 
for  processing  and  storing  input  data  from  the  video 
d^tizer  and  controlling  the  position  of  the  mechanical 
stage,  (8)  an  interactive  terminal  to  enable  operator 
control  of  the  system,  (9)  a  joystick  for  manual  control 
of  the  frame  position,  and  (10)  theXYZ  drive  ccmtndter 
to  move  the  holographic  frame  in  the  laser  beam.  Tabte 
I  presents  a  list  of  ^e  equipment  which  comprises  the 
system. 

The  system  is  designed  to  perform  the  following 
functions:  (1)  determine  the  position  of  a  particle  wi^ 
respect  to  a  3-D  frame  reference,  (2)  determine  the 
width,  height,  cross-sectional  area,  and  perimeter  of  a 
particle  using  a  scanning  algorithm,  (3)  determine  the 

TiMik  Mn«www»kimiaiii»iiiri^  -im  i  i 


Item 
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4 
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3 

Video  disitiMr 

Colorado  video 

270A 

6 
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Panaaonic 

WV6300 

7 

Micraeaaaputer,  Z-80, 
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Cromemoo 

Z-2D 

8 

CRTtafadaal 

Sorae 

iqi20 

8 

JoyalWi 

Croawnoo 

JB-IW 

10 

A  YA  QfIfV  QMnffOlHv 

Aantacli 

BC-2 

cross-sectional  area,  pmrimeter,  and  ed^  oowdinates 
of  the  p^cle  using  a  tracking  algcnitto,  (4)  enable 
particle  identifrcation  fixmi  frame  to  frame  1^  duplay- 
ing  the  particle  position  and  edge  points  (dupe)  fnmi 
the  previous  holographic  frame  on  the  current  frame, 
and  (5)  calculate  the  velocity,  average  diameter,  and 
particle  density  from  the  recmded  data  for  ea^  par¬ 
ticle. 

IV.  Teehnieal  Diaeuaaion 

Successive  holograms  of  settling  particles  were  re¬ 
corded  in  situ  on  negative  holographic  transparency 
film  using  the  particle  velocimeter  described  in  Ref.  6 
with  4X  magnification.  The  film  was  processed  using 
standard  photographic  techniques  and  mounted  on 
individual  frames.  Each  holographic  frame  was  posi¬ 
tioned  in  a  laser  beam  to  reconstruct  a  2-D  profile  of 
each  particle  in  the  direct  transmission  mode  with  an 
additional  magniflcation  of  25-40,  depending  on  the 
particle  size.  With  this  particular  setup  we  have  mea¬ 
sured  particles  down  to  10-/tm  diam  (see  Ref.  3).  In¬ 
creased  magniflcation  onto  the  holographic  fllm  is  re¬ 
quired  to  measure  smaller  particles. 

The  image  of  each  particle  was  successively  focused 
on  the  rear  projection  screen  by  moving  the  htdographic 
frame  in  the  Z  axis.  The  projected  scene  was  viewed 
with  a  TV  camera  and  converted  into  a  digitized  video 
signal.  The  holographic  frame  was  then  moved  in  the 
X  and  Faxes  to  bring  the  focused  image  of  each  particle 
into  the  center  of  the  TV  monitor  screen.  The  position 
of  the  frame  and  the  size,  shape,  and  orientation  of  ea^ 
particle  were  analyzed  and  recorded  by  the  computer. 
As  images  of  the  same  particles  in  subsequent  frames 
were  cataloged,  a  matrix  of  settling  velocity  and  lateral 
moticm  characteristics  was  generated  for  each  particle. 
The  statistical  distribution  of  settling  velocity  with  re¬ 
spect  to  particle  size,  shape,  and  orientation  provided 
an  accurate  determination  of  particle  density  and  ulti¬ 
mately  permitted  a  general  particle  classification  (e.g., 
organic,  mineral,  heavy  mineral),  based  largely  on 
density. 

A.  Particle  Position  and  Velocity 

Particle  settling  velocity  using  the  system  described 
in  Fig.  1  was  determined  by  measuring  the  position  of 
each  particle  in  3-D  space  at  two  or  more  different 
sampletimes.  A  tyincalpartide  configuration  is  shown 
inFig.2.  Prior  to  cataloging  any  particle  images  frxnn 
a  frame,  an  accurate  and  repeatable  referenoe  is  estab¬ 
lished  fw  a  point  common  to  aU  frames.  Tter^iNenoe 
VoisA  in  a  frame  is  centered  on  the  TV  monitor  by 
moving  the  frame  in  the  X,  F,  and  Z  axes  undw  com¬ 
puter  control  with  a  joystick.  The  image  of  eadi  par- 
tide  in  the  frune  is  tten  subsequmiti^  moved  to  the 
center  of  the  monitor  with  the  jctystick  and  its  positkm 
recordsd.  Sbe  and  orientation  parameten  an  deter- 
minsd  using  the  paitide  soann^  mad  sdge-tradi^ 
algorithms  described  below.  Onoereoordsd,tepar- 
tide  data  froBB  a  inevioui  fkaiae  can  be  nbfeved  and  tbs 
XYZ  stage  drivMi  to  that  poaitiQn  undir  eonmuftsr 
oontrol  to  enhance  looation  and  Manlttteaittoa  Of  that 
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opedaBMihip— ttdiiMqui>*i>lMchpfuc— wftiBvidM) 
ftamerfdata  to  obtain  adycooidiMtMafwiycoiiiplw 
•hapaa.  Howavw,  dua  nathod  irauld  ba  tor  rnoia 
ccnni^  and  tuna-conaumbu  than  naoaaaary  for  da- 
tarmioation  of  araa  and  pariaiatar  in  moat  of  our  ap- 
pUcationa.  Abo,  it  ia  rwa  for  nana  timn  ona  partida 
image  to  ba  found  in  a  given  focal  plana. 

Tliaiafora,  wa  have  choaan  an  imaga  acanning  tacfa- 
niqua  which  uaaa  a  aimpiified  edge  dataction  proceaa 
that  b  baaed  on  a  compariaon  of  the  piael  intanaity  with 
a  thraahold  adiuatod  for  the  avwaga  mtanaity  batwaan 
the  part^  and  backgroimd  to  detfnaiaa  pnrticb  area 
and  parimatar.  Pixab  at  each  tranaitton  acroaa  the 
thraahold  are  accumulated  aa  the  parimatar  of  the 
partida,  whib  pixab  above  the  thraahold  are  accmnu- 
Utad  aa  ita  araa. 


3.  Scamiiiixabaritim  pnoMMpmli  above  thfwhold  for  ifM 
and  parinMter. 


partida  on  a  aubaaquantfnma.  Whoi  the  entry  data 
for  all  framaa  haa  bean  completed,  the  dbtanca  aadi 
parti^  haa  traveled  batwaoi  aucc^ve  framaa  b  da- 
tarminad  and  the  individual  valodtiaa  cakubtad. 

B.  ParUda  Araa  and  Parimaiar 

DatarminatHm  of  the  projactad  crom-aactional  araa 
of  a  partida  can  ba  derived  from  the  raoonatructad 
imaga  by  counting  the  number  of  pbwb  inaida  the  par- 
tidaboundary.  toadBBilarnwnnar,tbaparinidarof 
the  partida  can  ba  derived  tonn  tba  munbm  of  phmb 
dong  the  edge.  Onaofthamajorprabbmaidthanoby 
imaga  b  datactioa  of  the  pi»da  iHUch  lapiaoMit  the 
edge.  Yaldmovdty  dbcuaiad  an  edge  datoction  tech- 
ni^^  bated  on  maximidng  the  likelihood  ratio  with 
a  ain^pb  afa^^pam  region-growing  algorithm.  TUa 
lalio  b  a  coi^pivbott  of  the  intanaitba  of  M^bor^ 
piiab  within  the  aama  ohjoct  to  naighhoring  piaob 
wtthfaitwodiftonntohiaeta.  ThabooMlaryaacbion 
b  braad  on  oomparibt.tha  ratio  a  pwdoflnad 
thmhold  mmI  aocmmdating  thoia  ndgMMdnoda  which 
aw  eoMMond  la  ba  within  a  dngb  ohbct  into  ana 
eonlipioaB  raglon.  HdCaoMd  AggNmadhavadtaal- 


C.  Partida  Area  Scan  Tadmiqiia 

The  partida  acanning  algorithm  davdopad  for  thb 
ayatam  b  baaed  on  aaimplifiad  aingb-pam  edge  detec¬ 
tion  proceaa.  The  aurfaca  araa  arid  a^  parimatm  of 
each  particb  are  datarminad  by  con^Muring  the  ddao 
intenaity  of  each  pixel  to  a  adactad  thra^old  vdua. 
The  technique  conabta  of  prooaaaing  a  a«iea  of  vnticd 
acana  which  horixontally  cover  the  image  aa  depicted  in 
Fig.  3.  The  horizontd  (equatorid)  and  vertkd  (polar) 
axea  of  the  p^de  are  determined  by  maimaliy  ad- 
juati^  a  verticdly  acanning  curaor  to  the  ed|^  d  the 
particle  when  the  ^orithm  b  firat  executed  for  each 
measurement  Thb  limits  the  ragitm  under  inveatiga- 
tion  to  the  immediate  environment  of  the  partida  utd 
eliminates  the  generally  out-of-focus  neighboring  par¬ 
ticles  from  consideration. 

Beginning  tm  the  left  side,  the  partide  b  acaimed 
from  ^  to  bottom  in  aynchronbm  with  the  TV  raster. 
The  video  intensity  (0-265)  of  each  pixel  b  compared 
to  a  threshold  vdue  cwresponding  to  the  intensity  at 
the  parUde  edge.  If  the  pixd  intensity  b  greater  than 
the  threshold,  ^  pM  area  b  assigned  a  unit  value  and 
the  column  area  b  incrMnantad.  Otherwise,  its  vdue 
bsero.  If  the  pixel  area  vdue  bdiffnent  from  that 
the  last  pixel,  the  horiaootd  perimeter  b  incranented. 
When  the  vertkd  scan  b  com|date,  the  difbrence  be¬ 
tween  areas  of  currant  and  previous  columns  b  added 
to  the  vertkd  parimatar.  The  scan  b  then  inaamanted 
to  the  r^t  and  the  procem  repeated  untfl  the  entire 
partide  has  bean  scanned.  The  aelectad  scan  consists 
ofNriowsvartiidlyaadAleColunHishoiinabdy.  The 
totd  araa  of  the  pwtkb  b  the  sum  of  all  colunui  areas 
scaled  mahown  hi  Bq.(l).  Sfawethaahapeofeachpiid 
b  not  square,  the  vertkd  coatyOMot  nMHt  ba  acabd  to 
provide  an  accurate  maaBuiemant  of  the  ana  and  pe- 
rbaetar.  Orsha|dy» 
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(2).  If  VpdJ)  ia  not  equal  to  VpiiJ  4-  1),  the  incio- 
OMBtal  horiu^  coaponont  HpiiJ)  •  1.  otharwiaa 
HpdJ) ■ ft 
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The  mqior  error  in  thia  mothod  la  hi^ily  dopMidaot  on 
the  orientation  of  angr  edge  aafaent  with  roq^  to  the 
acaiming  exit.  Ahhciu^  atraifht  edfoo  orthofonal  to 
the  acan  will  yidd  good  niooooranmita  d  perhnater, 
thoaa  aanM  Btiaight  edgao  oriented  diagonally  to  tiia 
acan  axia  can  depade  the  perimeter  maaauremant  fay 
aamuGhaiaO%(l-SQirr(l/2)].  Therafore.  a  different 
methodofdateraiiniiigperiinateririiidiwfflbeiiiaan- 
aitiva to partida orientation wae required.  Aloo,vari- 
ationa  in  the  piid  intmaity  around  the  partide  perim¬ 
eter  with  raapact  to  the  throahold  level  incraaaa  mea- 
aumnantunoartainty.  A  dyoandc  throahold  iriiich  can 
adapttovaryhigpartidea^bockgroundvideointeB- 
aitiaaerffliamrovethiaconditien.  Suchatachniquaia 
implanMaalad  far  an  edga-traddng  routine  deocribad 
fadow  that  pawridaa  data  for  detennination  of  partida 


Boundary  definition  uring  direction  and  curvature 
chaina  haa  bean  daacribed  by  Bcelaa  H  at.*  Thaoa 
tadudquaa  foaeially  deacxfbe  an  fanage  boundary  in 
tanaa  of  a  1-D  Hat  of  Midw  atarting  flrmn  an  erigfaa  on 
thahnaai  adpi  and  inrrwaanHiig  unifonnly  around  the 


dataaiaiaiiif  ahipa  fay  dafta^  critical  poteta  anMBd 
thabaaWhvy.'*  Thaaopointafaidudemoontieuitiaa 


if -AN 
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Fif.  4.  Tnckinc  •Itotitlim. 

in  curvature,  pointa  of  inflection,  curvature  maiima, 
interaactiona,  and  pointa  of  tangency.  Fourier  trana- 
fonna  have  boon  umd  aa  digital  flHm  to  amooth  the 
digitiied  boundariea  of  idanar  objacta.  Pavlidia  dia- 
cuaaaa  aeveral  algorithna  uaing  Fouriw  tranefomn 
coeffkianta  to  de&ie  detailed  ahape  characteriatica." 
Teague  haa  developad  aimilar  ahape  analyaia  tadmiquea 
uaing  Zmiike  Brnmanta.^^ 

Wa  have  ehoaan  a  boundary  diain  tedmique  to  de- 
tormina  partida  ahape  and  orientation  which  ia  aimilar 
to  the  mdhod  daacribed  by  Eodea.*  The  adge-traddng 
algorithm,  daacribed  bd«m,determinea  the  boemdary 
between  the  particle  and  biwkground  by  examining  a 
nnall  legiw  of  pixela  along  the  partida  edge.  Edge 
detactMm  ia  baa^  on  the  intenaity  gradient  between 
neighbOTing  pixela  rather  than  the  aimple  threabold 
technique  ua^  in  the  particle  ■canning  algorithm  de- 
acribad  above.  Aa  ea^  new  point  akmg  tte  partida 
edge  ia  determined,  the  tracing  pattorn  ia  thm  cm- 
tered  on  that  point  and  the  new  naighboriiood  ia  ex¬ 
plored.  Thia  inocam  oontinuea  until  the  partida  haa 
been  folly  dreunmavigatad.  Amoredetdiaddiacuaaiao 
ia  given  below. 

E.  Perticia  Edga-TiacMng  Tadnlqua 
The  partida  edge-tracking  algorithm  providaa  the 


The  edga-traddng  algorithm  ocana  acrom  the  partida 
edge  far  a  aarim  of  amtiB  dreuiar  pattama  wharo  i  ia  oiw 
of  the  Np  acaaa  around  the  part^  edge  mid^  ia  one  of 
thoM;pointainaacanaaaliowBinFig.4.  Tlwpointy 
in  the  acan  for  vdddi  the  plial  faBtanaity  V(tV)  oxoaoda 
thaprevienapfodintenrity  V(iJ  -  1)  Iqr  a  throahold 
value  V,  ia  daflnad  at  the  edge  for  the  eurreut  aonn 
pattami  and  ia  need  aa  the  centar  of  the  nart  acan  pat 
tarB(t4>l).  Thia procim centinum uatfl tho acaa ro- 
tumatethaBaighboriwedofthaatartiagpoint  The 
acmi  tadhia  R,,  number  of  polma  par  aean  N,t  and 
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Fig.  7.  VidaoaakpUtndaacnaatliaadtaofahi^oaiitiaat(tap)aBd 
low  coBttaat  (bottom)  partida. 


The  perimeter  of  each  particle  is  measured  by  accu¬ 
mulating  the  radial  distance  of  each  threshold  transition 
around  the  particle  edge  until  the  circunmavigatkm 
the  particle  is  complete.  Since  the  shape  ttfeadt  pixel 
is  not  square,  both  vertical  and  horixontal  components 
must  be  scaled  accordingly,  by  Fy  and  fj, ,  reflectively, 
to  provide  an  accurate  measurement  ^  the  distance 
around  the  particle.  The  algorithm  for  t^  function 
is 

Ny 

porimotor  •  lum  |(#V  •  r,)*+  (f,  •  T,)*)'^,  (3) 

i  ■  1 

where  Ty  "  R,  sin. 

Tf^Rt  cos, 

R,  •  scan  r^us,  and 

Np  •  number  of  threshold  transition  points  along 
the  edge. 

The  area  of  each  particle  is  determined  by  accumu¬ 
lating  the  tnqieioidal  area  defined  by  each  new  transi¬ 
tion  as  shown  in  Fig.  4.  The  distance  Xi  between  the 
particle  edge  and  a  vertical  reference  UiM  at  the  center 
of  the  screen  is  determined  for  eadi  trantitioo  point  and 
multiplied  by  the  vertical  dirtanoe  dy^  betwemi  two 
suoceesiveef^tea.  The  entire  area  of  the  pgffticle  is  de¬ 
termined  by  ■uccessively  adding  or  subtractirf  eadi 
incranentai  area  depending  on  the  vertical  diracthm 
until  the  particle  has  beni  oomidetdy  encircled  as  de¬ 
scribed  by 


I 

-  aum  (*<  +  *,>i)  •  dFi  •  Fyif*  L  W) 

<1*  a  1  I 


Pirtideidsntificntionflwm  frame  to  ttaase  is  achieved 
by  diHiloying  the  sda  |dsdi  of  a  sslMled  particle  OB  the 
nertmiM  and  eoaporing  h  with  the  shapes  of  the  new 
paitidee  for  eondetion. 


6,  reapectively.  These  errors  are  functions  of  the 
numbw  of  scan  points  (N.)  and  the  scan  radius  (fid.  as 
shown  in 


radial  error  ■  A.  •  |1  -  eim{2wlN,)\l2, 

Angular  error  ■  +/—2v/S,. 

The  mor  in  detensining  the  perimetm  of  a  smooth 
boundary  will  produce  a  sinall  positive  error  (04.,  2jB% 
for  N,  *  ao)  as  shown  in  Eq.  (6).  However,  the  algo¬ 
rithm  may  exhibit  a  large  negative  error  for  figmes  witii 
sharp  corners  becauee  of  tnmcation  produ^  by  the 
dwcrateacanprooeas.  The  difliinmioe  between  tte  ac¬ 
tual  and  measured  area  for  large,  smooth  particles 
should  be  small  but  may  become  negative  for  particles 
with  sharp  corners,  Bq.  (7): 

parimatar  arrar  ■  +Np  •  A.(l  *■  ooa2r/A(.),  (S) 

araa  arrar  -  +/-{Np  •  A,  *  aia2r/N|)/l  (7) 

Area  and  pwimeter  measuremNits  will  also  vary  as 
a  function  of  the  video  intensity  threshold  vtiue  for  the 
scan  mode  and  as  a  function  of  die  video  intensity  gra¬ 
dient  tlireohold  value  for  the  tracking  mode.  Thevideo 
amiditode  acroee  the  edge  of  an  image  with  vai^ng 
ooatneteiesbownkiFlg.7.  Urn  sffocto  due  to  variation 
of  the  threshold  vahiea  on  meaauiementa  of  a  Ugh- 
contraet,  white-on-blaA  square  (2  cm/side)  are  sum¬ 
marised  in  Table  D.  In  generating  them  data  the  video 
threshold  iatomdtyfor&aoanBMide  was  set  at  valuee 
from  23%  to  74%  of  the  mtensity  diflersoce  between  the 
square  and  the  baAground.wlieieae  the  Inteneitygra- 
dient  thrediold  (intensity  dlffersnce/pixti)  foe  the 
tracking  algorithm  was  vaM  from  M  to  Sl%  of  the 
total  intensity  dUiermme  between  dm  square  end  the 
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Fig.  8.  Photograph  irf' video  image  of  kolognuB  of  mica  flak*.  The 
video  amplitude  eech  inml  along  the  central  white  line  ie  shown 
by  the  relative  scale  on  the  right-hand  aide  of  the  screen.  Thelength 
and  width  of  this  particle  were  41  and  19  pm,  reapectively. 


background.  A  negative  threahold  gradient  (i.e.,  a 
bright-to-dark  transition)  was  used  to  determine  the 
edge  in  this  particular  application. 

The  range  of  threshold  effects  is  smaller  for  the 
tracking  mode  than  the  scanning  mode  becaiMe  the 
maximum  video  intensity  gradient,  generally  found 
about  halfway  betwemi  the  intensities  ^  the  ima^  and 
the  backipround  (see  Fig.  7),  cannot  be  exceeded  without 
the  algorithm  losing  track  of  the  particle  edge.  The 
threahold  intensity  gradient  value  for  the  tracking 
program  is  generally  selected  to  be  smaller  than  the 
smidlest  normal  edge  gradient  (bottom  curve  in  Fig.  7) 
found  along  the  parUde  edge.  This  insuies  that  diiique 
approaches  to  the  edge  do  not  reduce  the  measured  in¬ 
tensity  gradient  below  that  of  the  threshold.  This  re¬ 
duced  threshold  value  causes  the  tracking  program  for 
the  high-cmitrast  part  of  the  particle  (upper  curve.  Fig. 
7)  to,  in  effect,  underestimate  the  sixe  of  the  particle  fw 
a  white-mi-blMk  image  and  overestimate  tte  sise  of  a 
particle  for  a  Uack-on-white  image,  especially  for  par- 
tides  of  varying  contrast  along  the  edges.  Forthewfaite 
square  on  l^k  background  image  addressed  in  TaUe 
n,  this  resulted  in  paitide  sixe  estimates  that  were  % 
to  1  %  pixels  lower  per  edge  transition  (see  upper  curve 
in  Fig.  7),  depending  on  the  video  threshdd  intmisity 
gradirat  selected  or  a  1.2  to  3.6%  underestimate  of  the 
length  of  eadi  side.  This  underestimate  would  be  ofEwt 
tosmneextentbytheradialerroreffectonthepOTim- 
eter  ae  deecribed  by  Bq.  (6). 

V.  Fertemmiwe  lednedeii 

The  performanoe  of  the  scan  and  track  aigoritiuns 
was  eviJurted  using  white-on-black  geometrical  figures 
(aqimre,  equilideral  triamde,  and  cbde)  of  controlled 
eteandehapa.  The  fiiPMaare  typical  in  projected  eiae 
onto  the  TV  screen  of  magnified  holopaam  of  80- 
160-pmpartkloa.  Theydonotpuehthelteiteofreo- 
ohitlon  of  the  qrelem  einee  our  primary  objective  was 
to  teat  the  petformmce  of  varioue  image  anatmti  rigo- 
rithmeaa a fbnctko target eiee and aa^  Tbestaa- 


duds  for  comparison  were  accurate  to  within  ^  to  1 
IMxel/side  or  to  about  ^2.5%  for  the  smaU  figures  and 
±1.25%  for  the  large.  The  results  of  the  comparwHtt 
are  shown  in  Table  III.  Digitisation  error  ^>pears  to 
have  produced  a  greater  effect  m  accuracy  of  tte  mwll 
figures  with  complex  shapes  compared  with  the  larger 
ones.  A  digitization  error  of  1  pixel/edge  can  result  in 
errors  as  la^  as  ±3.6%  for  the  small  and  ±2.0%  for  the 
large  triangle  measurements.  However,  this  effect 
should  approach  zero  as  the  shapes  become  larger  and 
smoother. 

For  the  tracking  program,  most  errws  did  not  exceed 
the  accuracy  limits  of  the  standard  measurement. 
Digitization  errors,  together  with  the  uncertainties  in 
the  standards,  can  account  for  the  diffumces  between 
the  measured  and  standard  sizes;  incluskm  of  the  po¬ 
tential  of  radial,  comer-truncation,  and  thresludd  in- 
tensityAhreehidd  intensity  gradient  errors  insures  that 
the  maximum  tracking  error  measured  (3.96%)  falls 
within  potential  error  limits. 

For  the  scan  program,  all  the  areas  measured  were 
within  expected  uror  limits,  althou|d>  the  trian^  and 
circle  perimeter  errors  were  laigu.  Theperimeterer- 
rws  for  the  squares  were  quite  low  as  the  edgm  of  eudi 
square  were  aligned  with  the  vertical  and  horizontal 
axes  of  the  TV  camera.  Perimeters  in  the  scan  mode 
can  be  exaggerated  by  as  much  as  1  -  fw  lines 
skqnng  45*  to  the  vertical.  The  largest  manifimtation 
of  this  ssrratsd-edge  effect  was  found  for  tils  code,  with 
mors  than  a  ±14%  error. 

n  AgplaallantaMsiairami 

The  apitikation  of  these  ta^niques  was  dsnK»- 
strated  by  rsaxuBining  in  sttn  transmisBiaa  hokfiaam 
of  mica  flakes  ooUacted  during  the  seeing  expaftamat 
reported  in  M.  5.  Because  the  mica  flrdno  an  very 
thibi  and  tend  to  settle  at  soase  nide  to  the  food  plane, 
only  one  side  can  be  in  abrnp  focus  at  any  tbne. 
Therefore,  a  sHgiitiy  defbcueed  (lemer  contmri)  image 
uns  anriyaed  (Ffg.  8).  This  lowwr  contrast  iasage  is 

**  I*— TT  I*  virt  tt  m  1 1*  limit  timrm  aw 
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typical  ot  the  type  collected  from  irregularly  shaped 
partk^  oMnmmbr  found  in  seawater. 

A  comparison  betwemi  the  tracking  and  scanning 
algorithm  data  (TaMe  IV)  shows  that  the  areas  gener¬ 
al  are  within  <rfeach  other  for  these  mka  holo¬ 

grams.  However,  the  elliptical  image  of  the  third 
meaeuremMit  resulted  in  a  perimeter  that  was  11% 
higher  for  the  scanning  than  for  the  tracking  algorithm. 
The  Img  axis  of  this  image  was  diagonal  to  the  screen 
axes,  and  the  scanning  algmthm  overestimated  the 
perimeter  hy  serrating  the  diagonally  oriented  edges. 

VM.  temmary 

The  holographic  microvelocimeter  has  proven  to  be 
an  invalui^  to(d  in  studying  the  settling  dynamics  of 
partides  in  both  marine  and  laboratory  environments. 
The  reduction  of  data  from  this  instrument  has  been 
seasiautomated  using  a  microcomputer  to  control  two 
fdiases  of  the  analysis  and  provide  data  sUwage  and 
ratrievaL  Duringthefintphaseoftheanalydsaprs- 
dshmtnmslatinn  stage  is  controlled  to  bring  a  sequence 
of  holognphic  frames  into  register  for  the  3-D  mea- 
sioeBMnt  of  parride  displaoement  betwemi  frames. 
Control  of  this  stage  is  of  the  order  of  A0.0002  cm  w^ 
even  graatm  accurades  in  particle  settling  velocities 
br^  adiieved  by  increasing  the  hologram  exposure 
iatorvaL  ^  _ 

of  the  sko,  ah^Mi,  and  periaiatar  of  the  digitfand  video 
innai  of  suspended  psaddsswitkia  dm  hefogimn.  Two 
tecfosiqusa  have  been  devekmed  to  aseaeuse  the  mea 
aadparfanatorofanhaage.  Both  seedMide  have  been 
tested  with  etandsad  gemnetifo  as  woB  aa  hologfaphic 
and  are  aecurate  to  within  4-41%.  Tnefiret 
uriham  a  scanning  algorithm  based  on  absolute 
video  iatanei^  for  edge  detection.  AMMUghthealgo- 
ritfam  is  qdfo  foet  and  ie  ehle  to  seen  dm  area  and  pe¬ 
rimeter  of  a  partide  in  a  eingfo  paie,  aeentncy  ie  com- 

tie  amoD  oanot/ mil,  Nat/ If  jmenyitM 


promised  for  very  noisy  low-omtrast  partides  or  for 
pertides  with  long  diegi^  edges.  The  secoiMl  method 
uses  a  tracking  algorithm  bnmd  on  video  intensity 
gradient  for  edige  detection.  Thie  technique  not  mtly 
provides  accurate  area  and  perimeter  measurements  but 
also  produces  detailed  edge  coordinates  which  can  be 
used  f(«  complex  shape  idmtifkation.  Althou^this 
system  has  bm  appM  to  the  study  of  partide  settling 
dynamics,  the  techniques  develop^  could  also  be  a^ 
plied  to  other  pattern  recognition  and  image  analysis 
problems. 


This  effort  was  funded  under  Office  of  Naval  Re¬ 
search  contract  N00014-76-C-0639  to  the  Marine 
Sdenoe  Department,  University  of  South  Florida. 
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THE  HYDRAULIC  EQUIVALENCE  OT  MICA' 
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AasTCAcr  SwiKt  expcrinieatt  peifonned  on  tih-  to  Bf  nod -tiled  mica  flakes  with  a  holagnvMc  aai- 
oo-vdodaiMr  lavcalad  Aal  mica  it  the  bydnailic  eqaivaleal  of  qaartz  ffhema  havia(  iUmBWaw  a  bclar 
of  4  to  12  tiomt  tmaUer.  Mica  ia  the  veiy  fine  to  fiae-ttad  liica  hat  hcea  tndMoeally  aied  by  tadiaiwi 
ulogiMi  to  deiiaemr  teeat  of  depotitioe  or  aoodepatilioe  aad  |wmatial  wiaaoariat  of  fiaet.  aad  it  hen 
foaad  to  be  the  hydmulic  eqaivakal  of  tih-tized  paitiekt  bat  aol  of  day. 

Etperiamatt  abo  thowed  that  aiica  flaket  lead  to  lettle  at  atiealatioat  which  aee  aeither  peipendicular 
aorp^kl  to  the  gnvilatioaal  vector  aad  lead  to  geaerallyoiaiBtaia  their  orieataliaethraaghaut.  Eqaathau 
fiir  the  teotiag  of  a  diac  ia  Lenaaa  aad  olhen  (1974)  aad  that  developed  by  Komar  aad  Reiaien  (1978) 
am  ihaam  lo  be  aiathematically  thnilar  for  the  coane-tih  to  fiae-taod  laaget  of  diact  and  ate  adeqaale 
peedidort  of  lelllieg  ralet  of  mica  flakes.  A  comparitoo  of  the  hydnulic  eqiiivaleacy  of  quaitz  spheres  to 
coarm^*  through  fiae-saod-tiaed  aiica  flaket  is  presemed. 


INTRODUCTION 

Became  it  cleaves  into  flaky  paitkies,  sand- 
sized  mica  is  a  mineral  group  in  which  shape 
should  obviously  affect  settling  charaleristics 
and,  in  fact,  has  often  been  considered  to  be 
die  hydraulic  equivalent  of  silt  and  clay.  Nei- 
heaei  (1965)  noted  the  close  assockdon  be¬ 
tween  sand-sized  mica  and  the  clay  fraction  of 
Georgia  estuaries,  an  associatkxi  also  recog- 
toed  by  Bimenncblum  (1966)  off  Israel. 
and  others  (1968)  used  the  abundance  erf' mica 
in  the  125-  to  25()-|un  size  fraction  to  delineate 
areas  of  the  south^tren  United  States  conti¬ 
nental  margin  that  might  be  undergoing  win¬ 
nowing  or  deposition  of  fines,  a  process  which 
otherwise  wcwld  be  masked  by  die  dominance 
of  reworked  Pleistocene  sands.  A  similar  ap- 
proadi  was  used  by  Ad^obe  and  Stanly  (1972) 
on  the  Niger  Shelf.  Dc^  and  others  (1979) 
and  Park  t^  PiHicy  (IWl)  discuss  die  signif¬ 
icance  of  ndca  content  to  dqx)8itk)i^  and 
erariomri  systems  of  the  nlnde  continental 
margin  of  the  Eastern  United  Stales. 

Deapile  its  intuitive  widespread  oe  as  a  hy¬ 
draulic  aiudQg  of  finer  sized  se^neats,  no 
qwntitalive  evaluation  of  the  sedkaeaiologic 
(dwracteiistics  of  mica  has  yet  been  under¬ 
taken.  the  purpose  of  dds  paper  is  to  deter¬ 


MHHRnpi  fWWSB  f^WHy  fWIHS 


mine  as  far  as  possible  the  hydraulic  charac¬ 
teristics  of  mica  flakes  and  how  these 
characteristics  compare  with  those  predicted 
from  settling  equations  in  the  literature. 

APPROACH 

Our  approach  is  to  utilize  a  modified  holo- 
graphic  mictD-velociineter  developed  and  de¬ 
scribed  by  Carder  (1978)  and  Carder  aiul  Mey¬ 
ers  (1979).  Figure  I  shows  die  system  used. 
Transmission  holograms  are  simultaneously 
collected  along  the  vertical  and  horizontal  axes 
of  a  settling  cuvette.  A  reference  poim  in  the 
cuvkte  is  also  recorded  on  each  vertical  and 
horizontal  hologram  to  allow  translation  be¬ 
tween  the  two  axes.  Sequential  frames  at  ac¬ 
curately  timed  intervals  (from  0.5  to  3.0  sec¬ 
onds  dqiending  on  the  particle  size)  record  the 
setding  velocity  and  orientmion  as  well  as  par¬ 
ticle  size  and  shape.  The  images  are  recon¬ 
structed  by  placing  the  hologram  back  into  the 
laser  pndi  a^  refocusing  on  particles  in  any  of 
the  infinile  number  of  focal  planes  in  foe  set- 
ding  cuvette.  This  system  ofito  the  advantafe 
of  using  actual  sedinsentary  particles  in  foe  snut- 
10  clay-size  ranges,  thus  obviating  any  ariois 
font  may  be  inherent  fat  scale  mode^  qrs- 
lems. 

Mica  samples  were  choaen  ftum  sedhnent 
cores  ftom  foe  Eaaieni  Unitad  Smies  continen- 
ttisiopa  (Doyle  and  ofosR,  1979).  Sfo- to  fina- 
sand-siaed  mica  Itakaa  wore  added  to  a  arnaR 
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Flo.  1.— tMognpUc  micfD-velociiiiclerteiupiHedto 
meaawe  the  senling  velocity  of  the  mica  flaket. 


quantity  of  filtered  distilled  water  to  obtain  a 
sluny.  A  drop  of  slurry  was  picked  up  on  a 
fine  brush  and  the  drop  barely  touched  to  the 
top  of  the  miniacus  of  the  distiUed  filtered  water 
in  the  cuvette,  thereby  introducing  mica  fiakes 
into  the  measuring  apparatus.  Odier  methods 
of  sample  introduction,  including  use  of  an  eye 
dropper,  were  tried,  but  they  tended  to  set  up 
convection  within  the  cuvette. 

Resulting  settling  velocities  were  then  com¬ 
pared  with  die  dieoretical  setding  rates  gen¬ 
erated  from  the  equations  of  Lerman  and  oth¬ 
ers  (1974)  and  Komar  and  Reimers  (1978).  In 
die  calculations  of  setding  rates  a  density  for 
mka  of  2.9  g/cm^  was  used,  midway  in  the 
normal  range  for  biotite  and  muscovite,  Ex- 
aminadon  the  mica  by  standard  optical  tech¬ 
niques  showed  that  the  mica  fracdon  in  the  cores 
was  composed  of  muscovite  and  a  lesser  amoum 
of  biodte.  Finally,  we  compared  die  measured 
setding  rates  of  the  mica  with  those  for  quartz 
spheres. 

SHAtCS  AS  A  FACroa  OF  ORAIN  SETIUNG 

Sluqie  has  been  recogni^  as  an  impoftaat 
factor  in  the  arudysis  of  sediments  Iqr  hydraulic 
methods  (prindpally  setding  tube)  fo  over  l(X) 
years.  Gibbs  sod  odien  (1971),  Lerman  and 
odiers  (1974),  Komar  and  Reoiien  (1978),  nd 
Btezina  (1979)  have  raramarized  die  devd- 
opmmit  of  thinking  conoerning  tte  hydraulic 
importance  of  sha^  and  have  contributed  to 
the  fonanlation  of  equations  for  aerffing  ve¬ 
locity  that  take  grata  shape  hito  aocouat.  la  a 
aeries  of  recent  aiticles  Bdba  «d  Komr  (1981a 
and  b)  have  bcfun  to  examine  slnyie  tiltaw  of 
various  types  of  naturai  particles. 


Lerman  and  others  (1974)  modified  the 
equadoos  for  the  settli^  of  a  disc  cf  “no 
thickness*  in  the  Stoke’s  range,  developed  by 
Lamb  (1932),  Payne  and  PleU  (1960),  and 
Brenner  (19M),  iy  adding  a  term  for  disc 
thickness.  The  resulting  formulae  for  the  two 
major  orientations  of  fall  are; 


..  g(p.-p)qhi^ . 

^  - (edgewise)  (1) 

3.396  |i 

..  g(Pi  “  P)<W^  ^  ^  , 

U  - (broadside),  (2) 

5.1  p. 


where  U  =  setding  speed;  r  -  disc  radius;  h 
=  thickness;  qi,  =  h/r,  p  =  dyruunic  viscosity 
of  water;  p,  =  particle  density;  p  =  density  of 
water,  and  g  =  accekradon  due  to  gravity  (981 
cm/sec^). 

Komar  and  Reimers  (1978)  approached  the 
effect  of  shape  on  setding  velocity  by  scale 
modeling  wito  pebbles  in  a  glycerine  setding 
medium,  the  results  being  equivalent  to  quartz 
sand  and  silt  in  water.  They  found  experimen¬ 
tally  duu  the  Corey  Sluqie  Factor  (CSF)  intro¬ 
duced  by  Corey  (1949),  Malaika  (1949), 
McNown  and  hMaika  (1950),  and  McNown 
and  odiers  (1951)  gave  the  best  prediction  of 
shape  effects  for  the  p^les  they  used  in  their 
experiment.  Based  u^  their  experimental  re¬ 
sults,  Komar  and  Reimers  (1978)  developed 
an  empirical  formula  for  setdtog  velocity  in  the 
Stoke’s  range,  taking  into  account  shape  ef¬ 
fect; 


U  =» 


1  1 
)8pf(CFS) 


(P.  -  P)  gDj, 


(3) 


where  O,  -  (hO|D|)''^;  Di  »  intermediate  axis 
particle  diamirter;  D|  »  principal  axis  particle 
diameter,  CSF  «  h/  VD|Di;  h  «  small  axis 
preticle  diameter,  f(CSF)  >  0.946(CSF)~*'”* 
when  0.4  s  CSF  s  0.8;  and  KCSF)  -  2.18 
-  2.09  {CSF)  when  CSF  s  0.4. 

Ihey  fordier  found  that  they  ooukt  extend 
the  range  of  the  grain  sizes  fUlowing  Stoke's 
setding  to  RejmcMs  numbers  of  iqi  to  0.10  or 
for  pteticles  ^  up  to  too  pm,  and  dtey  were 
able  to  empirically  extend  their  dam  to  cover 
gram  sizes  up  tnrougn  pennies,  imr  pamcies 
like  ntica  where  h  is  sflMU  lelntive  to  Dk  and 
D,.  CSF  is  9^0.4  and  the  geneial  eqnadoa  (3) 
becomes: 
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U  »  (1/18  li  ) - 

'  (2.18 -2.09  (h//DP,)] 

•(p.-plgCliDA)*'*  (4) 

The  perticle  thidmeas  h  may  be  expressed 
as  a  rttio  widi  respect  to  D,  that  is,  D/1(X), 
D/SO,  D/2S  .  .  .  D/x,  and  dierefore  in  the 
ranges  in  which  we  are  working,  the  equations 
of  Leman  and  othen  (1974)  ^  equation  (3) 
«e  siinilar.  Fdr  examide,  let  0|  -  D|  =  D,  x 
»  6/h  »  2r/h,  and  ^  »  h/r.  Then  equation 
(4)  becomes: 


U 


1  1 


18ti 


^2. 18- 2.09 


•(p. 


-p)«(y) 


u 


^2. 18- 2.09 


.in 


~(p.-p)gD^  (3) 
3L_- _ p 

Stoke’s  Formula 


Our  measurements  of  h  for  mka  flakes  yielded 
a  range  (rf  h  *  0.007D  -  0.05  D.  For  x- 
h  »  1/0.01  >  100,  and  substituting  in  (S),  we 
obtain: 

U  -  0.021  ;^(p»-p)gry. 

18(1 

Lerman’s  formulae  are  of  the  goteral  fom 
^  8(P>-P)V^ 

C|1 

wherec  >  3.396  to  S.l.  Substituting  as  before 
weobtifai 


U 


g(p.  -  p)  O^ 


2c|Ui 

Msddplyiiv  by  18/18  yields 

(Pb-P)8D*- 

2cx  I8»i 


Stoke’i 


sronmila 


For  Lerman’s  formula  to  be  equii^eat  to  Ko¬ 


mar  and  Reimeis  at  x  =  100  we  can  set 


J18 

2cx 


==0.021 


18 

^  ~  2(100X0.021) 


c  =  4.28. 


This  value  lies  between  Lerman’s  two  values 
of  3.3%  and  S.l,  and  therefore  values  derived 
ftom  the  two  sets  of  finmulae  will  be  close  to 
each  other  over  the  ranges  we  are  examining. 

Assuming  Dj  =  D,  =  D  (where  D  is  an  av¬ 
erage  diam^r)  is  an  ideali^  case  which  all 
mica  flakes  do  not  meet.  Any  deviation  from 
D|  =  D|  will  obviously  cause  some  deviation 
in  the  (^ulation  of  p.  When  Di  =  D„  the 
parameter  x  =  D/h  beiromes  the  shape  param¬ 
eter  direcdy  relat^  to  the  C!orey  shape  factor. 

TUB  HYDRAULIC  EQUIVALENCY  OF  MCA 

In  our  experiments  mka  flakes  seldom  set¬ 
tled  eitho'  exacdy  broadside  or  edgewise.  They 
usually  assumed  an  orienlation  betiimen  the  two 
extremes,  neither  perpendicular  nor  parallel  to 
the  gravitational  vector  whkh  they  maintained 
throughout.  Stringham  et  al.  (1969)  also  showed 
that  particles  setding  at  low  Reynolds  numbers 
maintain  their  initial  orientation  which  need  not 
be  perpendicular  to  the  setding  direction. 

Figim  2  shows  the  experimentally  deter¬ 
mined  setding  velocities  of  mka  flakes  ploned 
on  the  family  of  curves  of  Lerman  and  others 
(1974)  and  Konur  and  Reimers  (1978)  for 
dikknesses  of  D/200.  D/SO,  D/20,  and  D/ 
10.  Komar  and  Reimer’s  (1978)  equation  be¬ 
gins  to  diverge  from  those  of  Lerman  (1974) 
only  at  thicknesses  of  D/20  and  greater.  The 
thkknesss/D  ratios  of  32  mka  grains  from  a 
split  of  the  slope  sediment  used  to  die  setditig 
velocity  experiment  were  determined  by  scat- 
lermg  some  quartz  beadi  sand  on  an  SEM  stub 
and  then  sprinkling  the  mka-rkh  slope  sample 
over  it.  Some  of  die  idea  flakes  landed  on  edge 
and  we  were  dms  able  to  measure  thkkness- 
to-D  ratios.  Thicknesses  ranged  ftom  0.7  per¬ 
cent  of  D  to  5  percent  of  D,  with  the  awnage 
at  2  percent  and  a  standard  deviation  I.IS 
peront.  Figure  2  also  shows  dwt  most  of  the 
velodties  of  mka  fldees  that  we  measured  fidi 
wMda  the  envehqie  ctesaed  by  the  theoretical 
curves.  Scatter  at  the  coarser  end  is  probtibiy 
moody  the  result  of  variation  in  dtickness. 
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KOMAR  AND  RIEMERS 


LERMAN  EDGE 
DATA  ROMTS  ■ 


DIAMETER  IN  MICROMETERS 

no.  2. — Letiuai  and  othen  (1974)  equation  for  the  settling  of  discs  both  broadside  and  edgewise  plotted  along 
with  that  of  Komw  and  Rcimers  (1978)  for  thickness  of  D/XO,  D/SO,  D/20,  and  D/IO.  Note  that  the  equations 
yield  values  which  «c  very  close  to  eaA  other  at  thicknesses  of  D/SO  or  less,  thicknesses  most  prevalent  for  mica 
fUkes.  X's  aic  values  of  mica  flakes  we  measured. 


qnmnetfy  (for  example,  1^  and  settling 
orientatkin,  because  some  paitkies  in  this  size 
nu^  are  beginning  to  exc^  the  Stoke’s  low 
Ri^iBolds  nwnber  constraint  on  the  State’s 
equation. 

At  the  ooniacr  end  of  the  scale,  particles  fall 
doaer  to  the  equations  of  Lerman  and  others 
(1974)  than  to  that  of  Komar  and  Reaneia 
(1978).  Hfure  3  shows  the  diameter  of  spher- 
feal  mnutz  particles  of  density  of  2.83  wlticfa 
amdi  aetile  at  rsiBS  of  the  variously  steed  mica 
flakee  datemitMd  ftom  foe  equations  of  Ler¬ 


man  and  others  (1974)  for  h  =  D/200,  D/SO, 
and  D/20.  Our  experiments  diow  that  mica 
ranging  in  diameter  from  62  |un  to  2S0|un, 
with  a  thickness  of  no  more  tiuui  S  percent  of 
die  sieve  diameter,  is  the  hydraulic  equivalent 
of  S-|un  to  82-|un  quartz  sphnes.  Very  fine, 
sand-sized  mica,  then,  is  the  l^draulk  equiv¬ 
alent  of  silt-steed  quartz  spheres  whereas  fine- 
sand-steed  mica  is  die  equiWleat  of  sOt-  to  v^ 
fine  sand-rized  spheres.  Coarae-sUt-siied  nrica 
flakes  are  the  eqidvalent  of  fine,  and  very  fine 
enartz  silts. 
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Fki.  3. — ^The  hydraulic  equivalency  of  mica  to  quaitz  spheres.  Mica  sizes  are  coaise  slit  to  fine  sand,  and  thick¬ 
nesses  are  D/200,  D/30,  and  D/20. 


CONCLUSIONS 

1)  The  equations  for  determining  settling 
velwity  of  a  disc  presented  by  Lmman  and 
others  (1974)  and  dwt  developed  by  Komar  and 
Reimers  (1978)  for  particles  where  the  Corey 
Shape  Factor  £0.4  are  similar  for  thicknesses 
j,/  of  D/50  or  less. 

^  2)  Expoiments  on  fine-sand-  to  coarse-silt- 

Y  sized  mica  flakes  using  a  hologiaphic  micio- 
velocimeter  showed  tlM  grains  most  often  do 
not  orient  themselves  perp^icular  to  the  flow 
field  as  they  settle  but  that  they  tend  to  settle 
at  some  orientation  between  brxMdside  and 
edgewise. 


3)  Mica  in  the  coarse-silt,  very  fine  sand, 
and  fine-sand  sizes  is  the  hydraulic  equivalent 
of  silt-  and  very  fine  sand-sized  quartz  spheres 
according  to  the  relationships  shown  in  Figures 
2  and  3. 

Contrary  to  previous  assumptions,  mica  in 
the  conrse-sih  to  fine-sand  sizes  is  not  the  hy¬ 
draulic  equivalent  of  clay. 
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